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Abstract 

We present experimental measurements of the cholesteric-smectic phase transition of suspensions 
of charged semiflexible rods as a function of rod flexibihty and surface charge. The rod particles 



'^ ■ consist of the bacteriophage M13 and closely related mutants, which are structurally identical to 

this virus, but vary either in contour length and therefore ratio of persistence length to contour 
length, or vary in surface charge. Surface charge is altered in two ways; by changing solution pH 
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f~^ I compared with existing theoretical predictions for the phase behavior of flexible rods and charged 



a single charged amino acid for a neutral one per viral coat protein. Phase diagrams are measured 
as a function of particle length, particle charge and ionic strength. The experimental results are 



C^ ■ rods. In contrast to the isotropic-cholesteric transition, where theory and experiment agree at 
high ionic strength, the nematic-smectic transition exhibits complex charge and ionic strength 

^ dependence significantly different from predicted phase behavior. Possible explanations for these 

o' 

O ■ unexpected results are discussed. 

>: 

X ■ PACS numbers: 64.70.Md, 61.30.St 



* current address: Department of Materials Science and Engineering, University of Illinois at Urbana Cham- 
paign, Urbana, Illinois 61801 



I. INTRODUCTION 

In a suspension of hard or charged rods, purely repulsive entropic interactions are suffi- 
cient to induce liquid crystal ordering. Theoretically, hard rods exhibit isotropic, nematic, 
smectic and columnar liquid crystal phases with increasing concentration [l|, |2, |3|. Unfor- 
tunately, production of hard, rigid, monodisperse rods is very difficult. Rigid and flexible 
polyelectrolyte rods, however, are abundant, especially in biological systems, which by nature 
lend themselves to mass-production. In this paper we will study the influence of flexibility 
and electrostatic interactions on the formation of a smectic phase from a nematic phase 
using suspensions of charged, semiflexible fd and M13 virus rods. Viruses, such as fd, M13, 
and Tobacco Mosaic Virus are a unique choice for use in studying liquid crystal phase be- 
havior in that they are biologically produced to be monodisperse and are easily modified 
by genetic engineering and post-expression chemical modification. These virus particles and 
/9-FeOOH rods are, to our knowledge, the only colloidal systems known to exhibit the pre- 
dicted hard-rod phase progression from isotropic (I) to nematic (N) or cholesteric and then 
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to smectic (S) phases with increasing rod concentration |^, |5|,|6[. Even though qualitative 
theories have been developed to describe either the effects of electrostatics or the effects of 
flexibUitv on the ne.at.c-s.ectic ,N-S) pK.e transition of hard rods 003, tKey Kave 
yet to be thoroughly tested experimentally. Near the N-S transition, the particles are at 
very high concentrations, and as we will show, dilute-limit approximations of interparticle 
interactions, which are appropriate at the isotropic-nematic transition, cannot be used. By 
measuring the N-S transition of charged and/or flexible rods we learn about both the influ- 
ence of these parameters on smectic phase formation and the interactions between these rods 
in concentrated suspensions. Additionally, our results add insight into the ordering of other 
important rodlike polyelectrolytes such as DNA, which often appears in high concentrations 
under physiological conditions and exhibits cholesteric and columnar, but not the smectic. 
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In this paper we test the limits of current theoretical predictions for the nematic-smectic 
phase transition in three ways. First, we measure the phase transition for semiflexible ffi- 
amentous virus of identical structure and varied length. By changing the rod length and 



leaving local particle structure constant, the persistence length P of the rods, defined as one 
half the Kuhn length, remains constant. Subsequently, the rod fiexibility, as defined by the 
ratio of persistence length to contour length L, or P/L, is altered. In our experiments the 
fiexibility of the particles remains within the semifiexible limit, meaning P r^ L. Altering 
the particle fiexibility within the semifiexible limit probes the competition between rigid and 
fiexible rod phase behavior. Second, we vary the ionic strength of the virus rod suspensions 
allowing us to probe the efficacy of theoretical approximations for incorporating electrostatic 
repulsion into hard-particle theories. Third, we measure the nematic-smectic phase transi- 
tion for filamentous virus of different charge. Altering the surface charge by two independent 
techniques, solution chemistry and surface chemistry, probes the importance of the details 
of the surface charge distribution in determining long range interparticle interactions. By 
varying these three independent variables, length, charge and solution ionic strength, we 
systematically examine how electrostatic interactions and flexibility experimentally effect 
the nematic-smectic phase boundary. 

The colloidal rods we use are the rodlike semifiexible bacteriophages fd and M13 which 
form isotropic (I), cholesteric (nematic) and smectic (S) phases in solution with increasing 
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virus concentration [a 112, 11^, 11^ ■ The free energy difference between the cholesteric and 
nematic (N) phases is small, and therefore it is appropriate to compare our results with 
predictions for the nematic phase |l5|. Furthermore, near the nematic-smectic transition, 
the cholesteric unwinds into a nematic phase [lC)j. M13 and fd are composed of 2700 major 
coat proteins helicoidally wrapped about the single stranded viral DNA. They differ from 
one another by only one amino acid per major coat protein; the negatively charged aspartate 
(aspi2) in fd is substituted for the neutral asparagine (asni2) in M13 |l3]. They are thus 
ideal for use in studying the charge dependence of the virus rod phase transitions. Changes 
in the surface charge of the particles were also achieved by varying the pH of the solution 
|lq . Additionally, by varying the length of the M13 DNA we created M13 mutants which 
differ only in contour length. The M13 mutants have the same local structure, and thus we 
assume persistence length, as M13. These mutant M13 viruses were used to measure the 
flexibility dependence of the nematic-smectic phase transition. 



II. ELECTROSTATIC INTERACTIONS 

For colloidal rods, the total rod-rod interparticle interaction includes a combination of 
hard core repulsion and long ranged electrostatic repulsion. We present here two ways 
which have been previously proposed for incorporating electrostatic interactions into hard- 
rod theories for the nematic-smectic phase transition. The first originates from Onsager's 
calculation of an effective hard-core diameter (-Deff) which is larger than the bare diameter 
D. Dcs is calculated from the second virial coefficient of the free energy for charged rods in 
the isotropic phase [ij. Specifically, for hard, rigid, rodlike particles, the limit of stability 
of the isotropic phase against a nematic phase is given by the Onsager relation bci = 4 

n 

where b = ttL'^D/A and q is the isotropic number density of rods 1]. For charged particles, 
Onsager showed that the stability condition remains unchanged provided D is replaced with 
Des, thus besCi = 4, with bes = i^L'^Des/4:. Increasing ionic strength decreases D^s, and 
for highly charged colloids, like M13 and fd, D^s is nearly independent of surface charge 
due to the non-linear nature of the Poisson-Boltzmann equation, which leads to counterion 
condensation nea. the col.o.d snrfaee Q. In previous wo* QQ. the predict.on that ... 
is constant at the I-N phase boundary has been experimentally verified at high ionic strength 
(I> 60 mM, large L/Dcs) for our system of virus rods, as illustrated in Fig. ^ These results 
validated the mapping of the I-N transition of charged virus rods at high ionic strength onto 
a hard rod theory by using an effective hard diameter, Dgfr- At low ionic strength (I< 60 
mM), the prediction that feeffCj is constant did not hold because of the breakdown of the 
second virial approximation at small L/Des |l9| . Consequently, we expect that if D^s works 
to describe the N-S transition of virus rods it would most likely be at high ionic strengths. 
Since Dcs is valid only where the second virial approximation is valid, ie. in the isotropic 
phase, Stroobants et. al developed an approximate way to describe the electrostatic inter- 



actions in the nematic phase using this second virial approximation 2l|. They defined a 
nematic effective diameter D^, which is calculated from the isotropic effective diameter D^^: 
D,Nff = Deff[l + /^r7(/)M/)], where 

p(/) = l((sin0)) (1) 
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FIG. 1: (Color online) Isotropic-Nematic phase transition ^effQ plotted as a function of D^s for 
both M13 and fd suspensions in Tris-HCl buffer at pH 8.2. The original data for this figure was 
published previously 19]. The solid line is the hard-rod prediction for semiflexible rods with a 
persistence length of 2.2 //m|20|. For small values of Dgff (high ionic strength), the coexistence 
concentrations for the charged rods are effectively mapped to the hard-rod predictions. The ionic 
strength scale is for fd suspensions (M13 has a lower surface charge, thus D^s at the same ionic 
strength is slightly larger). 



and 



Vif) = -((-sm01og(sm0))) - (log(2) - l/2)p(/) 
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The average (...) is over the solid angle Q weighted by the nematic angular distribution 
function /(fi) with (p describing the angle between adjacent rods J22]. The parameter h = 
K,~^/Dcs, where k,^^ is the Debye screening length, characterizes the preference of charged 
rods for twisting. Crossed charged rods have a lower energy than parallel charged rods, 
and h correspondingly increases with increasing electrostatic interactions (decreasing ionic 
strength). This definition for the nematic effective diameter is accurate as long as the average 
angle between the rods and the nematic director yiff^ is much greater than D^^/ L [23]. In 
this limit, the second virial coefficient is still much larger than the higher virial coefficients, 
which can be neglected. Near the N-S transition the order parameter, as determined by 



x-ray measurements of magnetically aligned samples, is S = 0.94 



24| . Using an angular 



distribution function with this order parameter oi S = 0.94 we find that D^ = l.lGDes 
at 5mM ionic strength and D^g = l.lODcs at 150 mM ionic strength. This corresponds to 
^/(W) ~ QD^s/L for the largest value of D^^. 

Previously, we argued that D^^, which is independent of virus concentration, could de- 
scribe the electrostatic interactions of fd virus suspensions at the nematic-smectic transition 
J8|. However, as mentioned above, the use of Z^efr beyond the regime where the second virial 
coefficient quantitatively describes the system is not justified, and from our calculation of 



yiP^ we know the use of the second virial approximation is questionable. In this article, 
our expanded range of measurements of the N-S transition of virus suspensions as a function 
of ionic strength, virus length, and virus surface charge will demonstrate that electrostatic 
interactions at the nematic-smectic transition are much more complex than those predicted 
at the limit of the second virial coefficient. 

An alternative method for incorporating electrostatics into a hard-rod theory for the N-S 
transition was developed by Kramer and Herzfeld. They calculate an "avoidance diameter" 
-Da which minimizes the scaled particle expression for the free energy of charged parallel 
spherocylinders as a function of concentration 7]. With respect to ionic strength, D^ ex- 
hibits the same trend as D^s-, but unlike Z^off, -Da is inherently concentration dependent, 
decreasing with increasing rod concentration. D^. can never be greater than the actual rod 
separation, something which is not impossible with Deff- Furthermore, by using the scaled 
^article theory, third and higher virial coefficients are accounted for in an approximate way 



25l |26| . unlike Onsager's effective diameter. This makes the "avoidance diameter" more 
appropriate for incorporating electrostatic interactions at the nematic-smectic transition. 
One disadvantage with the free energy expression developed by Kramer and Herzfeld is that 
it does not reduce to Onsager's theory in the absence of electrostatic interactions. Never- 
theless, Kramer and Herzfeld's calculations do qualitatively reproduce previously published 
data for the N-S transition of fd virus [7]. However, the limited range of data previously 
available did not include some of the interesting features described in this theory, which we 
are now able to test. 
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III. MATERIALS AND METHODS 

Properties of the wild type {wf) virus fd and M13 inc 
diameter D = 6.6 nm, and persistence length P =2.2 /im ^^ 
same diameter as the wild type M13 and lengths of 1.2 /im, 0.64 /im, and 0.39 /im |l4| . 
Because the molecular weight of the virus is proportional to its length, the molecular weight 
of the M13 mutants is M = MyjtL/ L^t, with M^j = 1.64 x 10'' g/mol and L^t = 0.88 /xm, the 
molecular weight and length of wild type M13, respectively. Virus production is explained 
elsewhere J27|. Two of the length-mutants (0.64 /im and 0.39 /im) were grown using the 
phagemid method IJ, |27[, which produces bidisperse solutions of the phagemid and the 
1.2 /im helper phage. Sample polydispersity was checked using agarose gel electrophoresis 
on the intact virus, and on the viral DNA. Excepting the phagemid solutions which were 
by mass 1.2 /im helper phage, virus solutions were highly monodisperse as indicated by 



sharp electrophoresis bands. All of these virus suspensions form well defined smectic phases 

Q. 

All samples were dialyzed against a 20 mM Tris-HCl buffer at pH 8.2 or 20 mM Sodium 
Acetate buffer adjusted with Acetic Acid to pH 5.2. To vary ionic strength, NaCl was 
added to the buffering solution. The linear surface charge density of fd is approximately 
10 e~/nm f3.4± O.le^/coat protein) at pH 8.2 and 7 e~/nm (2.3± O.le^/coat protein) at 
pH 5.2 [18]. The M13 surface charge is 7e^/nm (2.4± O.le^/coat protein) at pH 8.2 and 
3.6 e~/nm (1.3± 0.1e~/coat protein) as determined by comparing the M13 composition and 
electrophoretic mobilities to those of fd |l9|. On the viral surface, fd has four negatively 
ionizable amino acids and one positively ionizable amino acid per coat protein. At neutral 
pH, the terminal amine contributes approximately +l/2e charge. The M13 surface has three 
negatively ionizable amino acids, one positively ionizable amino acid and the terminal amine 
(+l/2e) per coat protein. 

After dialysis, the virus suspensions were concentrated via ultracentrifugation at 200 OOOg 
and diluted to concentrations just above the N-S transition. Virus suspensions were then 
allowed to equilibrate to room temperature. Bulk separation of the nematic and smectic 
phases is not observed, perhaps because of the high viscosity of the suspensions near the 




FIG. 2: (a) Differential interference contrast microscopy image of nematic-smectic coexistence of 
fd virus suspensions. The smectic phase can be recognized by the ladder like structures. The virus 
rods are oriented perpendicular to the layers as illustrated. The uniform texture is the nematic 
phase, (b) Digitally enhanced image of (a). The scale bar is 10 /um. 

N-S transition. However, smectic or nematic domains can be observed using differential 
interference contrast microscopy in coexistence with predominantly nematic and smectic 
bulk phases, respectively. Typically coexistence is observed as ribbons of smectic phase 
reaching into a nematic region, as shown in Fig. |21 

The location of the N-S transition was determined by measuring the highest nematic 
volume fraction (0^) and the lowest smectic volume fraction {(p^) observed. The volume 
fraction 0"^ = c^v, where c^ is the number density, and v is the volume of a single rod 
7iLD^/A. A similar equation holds for (p^ . The concentration of the phases was measured 
by absorption spectroscopy with the optical density (A) of the virus being A269nm = 3.84 
for a path length of 1 cm. 

Since knowing the surface charge of the virus is critical to our analysis of the N-S tran- 
sition, we experimentally measured the pH of the virus solutions at concentrations in the 
nematic phase just below the N-S transition. We found that for an initial buffer solution at 
pH 8.2 (Tris-HCl buffer pKa= 8.2), the pH of the concentrated virus suspensions is slightly 
less than 8.2, but still well within the buffering pH range (pH=pKa±l). The surface charge 
does not change significantly over this range [l8[. At pH 5.2 (Acetic acid buffer pKa= 4.76) 
the measured pH of the virus suspensions near the N-S transition was slightly higher than 
5.2, with pH increasing slightly with decreasing ionic strength, most likely due to the rela- 
tively high concentration of virus counterions (50-100 niM) as compared to buffer ions (20 
mM). This shift further away from the pKa may influence the phase behavior by increasing 
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FIG. 3: Volume fraction at the nematic-smectic phase transition, (p , for multiple ionic strengths 
at pH 8.2 as a function of rod length L and flexibility L/P. On the right axis is the measured 
concentration in mass density p^ = (j) M/vNa, where Na is Avagadro's number. Legend for ionic 
strengths is as follows: o 5 mM, ▲ 10 mM, • 60 mM, v 110 mM, ■ 150 mM. With increasing 
ionic strength i;^'^ increases due to increasing electrostatic screening. Dashed lines are a guide to 
the eye at constant ionic strength. Within experimental accuracy the smectic phase transition is 
independent of flexibility within the range 0.18 < L/P < 0.54. 

the the viral surface charge. The imphcations of these measurements are discussed in the 
Results section. 



IV. RESULTS 



A. Flexibility and ionic strength dependence of the N-S transition 

Fig. El shows (f)^ as a function of the M13 mutant particle length, and therefore virus 
flexibility by L/P, for multiple ionic strengths. Focusing on how rod flexibility influences 
the phase transition, we observe that at each ionic strength the measured 0"^ is independent 
of virus length, within experimental accuracy, and thus independent of changing flexibility 
in the range of 0.18 < L/P < 0.54. The bidispersity of the 0.64 /im and 0.39 ^m rod 
suspensions does not seem to influence the phase boundary because these samples, which 
are 20% 1.2 /xm rods by mass exhibit the same phase behavior as samples which are 100% 
1.2 fira rods. 
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FIG. 4: Average values of <^^ (open) and 0"^ (solid) at the N-S transition as a function of ionic 
strength at pH 8.2. Average at each ionic strength is over the results for the four M13 length 
mutants. The solid line is ^^ taken from simulations by Kramer and Herzfeld |7| for the N-S 
transition of particles the same size as jd and with a renormalized surface charge of le~/7.1 A. 
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As the N-S phase transition is independent of rod flexibihty for these experiments, we 
averaged the results for ^"^ and 0^ from all particle lengths to study the ionic strength 
dependence of the phase transition. These averaged values for (^^ and 0^ are shown as 
a function of ionic strength in Fig. |3] We observe that with increasing ionic strength, 
and therefore a corresponding decrease in electrostatic interactions, the volume fraction 
of the phase transition increases until an ionic strength of approximately I = 100 mM, 
at which point the phase transition becomes independent of ionic strength. Whereas the 
increase in phase transition concentration with ionic strength has been observed previously 
for suspensions of /rf virus [s|, the plateau in 0"^ and 0^ at high ionic strengths is previously 
undocumented. 



B. Surface charge dependence of the N-S transition 

To determine the influence of virus surface charge on N-S phase transition, we measured 
phase behavior of jd and M13 at both pH 8.2 and pH 5.2. Fig. El presents the ionic 
strength and pH dependence of the N-S phase transition for jd (a) and M13 (b). Below 
100 mM, there is a strong pH dependence in the N-S transition, as shown in Fig. EKjb. 
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FIG. 5: Neniatic-sinectic phase transition volume fraction i;^ as function of ionic strength for 
suspensions of a) fd and b) M13 at pH 5.2 (sohd) and pH 8.2 (open). Figure c) shows M13 (pH 
8.2) and fd (pH 5.2) at 7 e/nm surface charge. Sohd hues highhght the ionic strength independence 
at high ionic strength. 

Suspensions at higher pH (higher surface charge) consistently enter a smectic phase at lower 
concentrations. Above about 100 mM, 0"^ is independent of ionic strength, as in Fig. |3 
The fd phase boundary at high ionic strength saturates around (p^^^ ~ 0.21, independent 
of pH, and the M13 phase boundary saturates around 4>^^^ = 0.24, also independent of pH. 
Even at these high ionic strengths, where the surface charge of the virus is well screened, 
the higher charged fd suspensions have a phase boundary at lower concentrations than the 
Ml 3 suspensions. 

In Fig. Et we compare the phase transition for M13 and fd suspensions when both viruses 
have the same surface charge of 7e~/nm. Because the rods have the same surface charge, 
we assume the rods differ only by the location of the charges (positive and negative) on the 
surface. At low ionic strength the phase behavior is similar, but fd suspensions consistently 
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enter the smectic phase at a shghtly higher concentration. We note that the small measured 
increase in pH at low ionic strength of the pH 5.2 viral solutions mentioned in the Materials 
and Methods section does not account for this difference. An increase in pH would lower, 
not raise, the pH 5.2 phase transition concentrations by increasing electrostatic interactions. 
At high ionic strength, the reverse is true; fdhas a lower phase transition concentration than 
M13 suspensions, as mentioned above. We believe the measured differences in (p^^^ between 
M13 and fd to be statistically significant, and will discuss this unexpected observation further 
in the following section. 

V. DISCUSSION 

A. Flexibility and ionic strength dependence of the N-S transition 

The nematic-smectic transition of flexible, hard rods has been studied both theoretically 
and computationally [3|, |9|, |28| . A small amount of flexibility is expected to drive the smectic 

Ehase to higher concentrations, from the predicted hard-rigid-rod concentration of (f)'^=0A7 
I, to approximately 0.75 < (p^ < 0.8 within the semiflexible limit J2|. Within the semiflex- 
ible limit {L/P ~ 1), however, (p^ is predicted to be essentially independent of flexibility 
J9|. This insensitivity of 0"^ to flexibility in the semiflexible limit is in agreement with the 
measurements presented in Fig. El We note that this result is in striking contrast to the sig- 
niflcant flexibility dependence measured at isotropic-nematic transition for this same system 
of semiflexible M13 mutants which we describe in a separate report |19l |. 

As our rods are charged, the ionic strength of the virus suspension plays an important 
role in determining the phase boundaries by screening electrostatic interactions. To compare 
our charged-flexible-rod results with current predictions for the N-S phase transition of 
hard (rigid or flexible) rods, we have to effectively account for the electrostatic interactions 
between our virus rods. Two methods for incorporating electrostatics into the N-S phase 

unn 

transition of hard rods [21, 122, 12^ were presented earlier in this paper. One way to do this is to 
graph 0fg, the measured effective volume fraction along the nematic-smectic transition, and 
compare it to the theoretical volume fraction, (p^^, for the N-S transition of hard, semiflexible 
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FIG. 6: Effective volume fraction volume fraction along the nematic-smectic phase transition 
'^cS ~ '?^'^(^eff)^/-^^ ^o^ multiple ionic strengths at pH 8.2 as a function of L. (p^ , the actual 
volume fraction at the N-S transition, is shown in Fig. El Legend for symbols is to the right of 
the figure. Dashed lines drawn are a guide to the eye and are at constant ionic strength. Because 
0fg strongly depends on ionic strength, we conclude that D^g does not describe the electrostatic 
interactions at high virus concentrations. 

rods la y. (pes is defined as c^7cL{D^fif/A = (p'^Dff^f / D^ and is shown in Fig El If the 
effect of electrostatics can be accounted for by replacing D with D^, as can be done at the 
isotropic- nematic transition at high ionic strength, we could predict that 0fg = 0fj^ = 0.75. 
In other words, if D^^ accurately models the interparticle electrostatic interactions, the 
effective volume fraction 0fg should be equivalent to the hard-flexible rod volume fraction 
and should be independent of ionic strength. Thus multiplying the measured values for (p^ 
shown in Fig. Oby (D^)^/D^ should result in the collapse of all the different ionic strength 
data. 

However, we find that <^fg depends quite strongly on ionic strength in Fig. IHl Previously, 
we observed 0^^ = 0.75 independent of ionic strength for suspensions of fd virus js|. The 
data in Fig. IHlis consistent with this value at high ionic strengths (60mM < I < 150 mM), but 
by including a larger range of ionic strengths as well as multiple particle lengths, we clearly 
observe an ionic strength dependence in 0^, with 0^ ranging from 2.5 to 0.5. Furthermore, 
the plateau in the phase transition concentration at high ionic strength is not captured by 
scaling (j)^ by (D^)^/D^, as shown in Fig. HI by the dashed curve for 0fg = 0.75. The 
large ionic strength dependence of 0^, and the measured ionic strength independence above 
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100 inM indicate that D^^ is inadequate for describing the electrostatic interactions at the 
N-S transition. This is in contrast to the I-N transition, where Des accurately incorporates 
the electrostatic interactions between virus rods at high ionic strength [ll^. Furthermore, 
because ^f^ > 1 at low ionic strength we conclude that D^ overestimates the electrostatic 
interactions. This is not surprising because D^ is based on the second virial approximation 
which, strictly speaking, is valid only for isotropic suspensions at low concentrations. Using 
Dcs to relate the phase behavior of charged rods to hard-rod predictions is inappropriate at 
the N-S transition, particularly at high ionic strengths. 

The plateau in (p^ at high ionic strength is indeed predicted for parallel, charged, rigid 
spherocylinders with a concentration dependent avoidance diameter D^^, as shown by the 
solid line in Fig 01 This avoidance model predicts that at high ionic strength (p^ saturates at 
'/'fat = 0.47, the theoretical value for hard, rigid spherocylinders [2, llf- Because our rods are 
semiflexible, we would correspondingly predict that 0f^^ would be equal to the theoretical 
value for hard- semiflexible rods, (p^^ = 0.75. Instead of this value, our measurements of the 
phase transition volume fraction as a function of ionic strength yield 0|^^ = 0.21 — 0.24, 
which is three times lower than predicted by hard flexible rod theories. This suggests that 
either the flexibility of the rods has lowered the phase transition from that of hard rods, in 
contradiction to both theories and simulations, or that the electrostatic interactions between 
the rods are not accurately represented by either Da or Dcs- We will return to this question 
at the end of the Discussion section. 

B. Surface charge dependence of the N-S transition 

The pH dependence of the phase transition visible in Fig. and the difference between 
Ml 3 and fd saturation concentrations present even when they share the same surface charge 
(Fig. Efc) also indicate that neither D^.^ nor Da are appropriate for describing the electrostatic 
interactions between the virus rods at the N-S transition. The non-linearity of the Poisson- 
Boltzmann equation predicts that for high linear charge density the long-range electrostatic 



potential between rods is insensitive to surface charge changes and thus pH changes 



fl 



2l| 



This is conflrmed at the isotropic-nematic transition, where the charge dependence is well 
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described by D^s and the pH dependence of the phase transition is very small [l9|. However, 
a strong pH dependence of the smectic phase transition at low ionic strengths for both fd (Fig. 
Et-) and M13 (Fig. Eb) suspensions is observed. The observed difference between M13 and 
fd saturation concentrations at high ionic strength and equal surface charge (Fig. Et) is also 
not expected from Poisson-Boltzmann theory. This high sensitivity of the N-S transition to 
changes in pH and surface charge configuration indicates that the charge independent nature 
predicted by both D^d and Da does not correctly characterize the electrostatic interactions 
at the concentrations of the N-S transition. 

One possible explanation for why the high ionic strength, and correspondingly high con- 
centration, phase behavior is sensitive to surface charge configuration is that the adjacent 
virus surfaces are separated by approximately one virus diameter (6.6 nm) [2^, which is 
on the order of the spacing between viral coat-proteins (1.6 nm), and the Debye screening 
length K = 3.0A/v^ = 9A. When the surface to surface distance is of the order of the 
Debye screening length, the continuous charge distribution approximation used in Poisson 
Boltzmann theory can no longer be used, as is done in the effective diameter calculations. 
Furthermore, it has been shown theoretically that discretization of the surface charges can 
change the predicted counterion condensation from that predicted by the non-linear Poisson 
Boltzmann equation J29|, |30|. Perhaps it is because we are in the regime where the sur- 
face charge configuration can no longer be neglected that we observe charge-configuration- 
dependent saturation of the nematic-smectic phase transition. Previous work by Lyubartsev 
et. al. has been done to simulate the electrostatic interactions between these virus rods in 
the presence of divalent ions using an approximate discrete charge configuration J3l|. We 
propose that theoretical models or simulations similar to those by Lyubartsev et. al of the 
electrostatic interactions of a dense, rod-like polyelectrolyte system which include the detail 
of the surface charge configuration including the location of positive and negative amino 
acids on the viral surface, may shed light on the experimental differences between M13 and 
fd nematic-smectic transitions at high ionic strength. 
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C. Origin of the ionic strength independence of the NS transition 

At high ionic strength we have measured an ionic strength independent N-S transition. 
This is similar to the phase behavior predicted by Kramer and Herzfeld, in which the phase 
transition volume fraction saturates at that predicted for the N-S transition of hard rods. 
Yet, the value for our measured (p^^^ is three times lower than the predicted N-S transition 
volume fraction for semi-flexible rods 0^ = 0.75. These observations point to a failure of 
theory to describe the role of electrostatics and/or flexibility on the N-S transition. 

It has been shown by Odijk that undulations of semifiexible rods in a hexagonal con- 
figuration are typically contained within a tube of a diameter larger than the bare rod 
diameter J32]. These undulations create an effective repulsion similar to Helfrich repulsion 
of membranes. The dominant interaction between charged flexible rods in a dense suspen- 
sion therefore depends on the relative size of the tube diameter and the electrostatic effective 

n 

diameter [2^. At low ionic strength the interparticle interactions would be dominated by 
electrostatics, and at high ionic strength the interparticle interactions would be dominated 
by steric interactions of the flexible rods. This simple argument qualitatively agrees with 
our observed phase behavior. If charged, flexible rod behavior is indeed dominated by steric 
interactions at high ionic strengths, we expect that the N-S phase transition would be at 
a value equal to that predicted by theory for hard-semiflexible rods. However, both theory 
and simulations for hard-flexible rods which incorporate this repulsion due to flexibility, pre- 
dict that flexibility destabilizes the N-S transition, subsequently increasing, not decreasing 
as measured, the transition concentration above that predicted for rigid rods jjjlsl, |3J]. It 
is possible that the role of flexibility is not accurately incorporated into the theories which 
predict an increase in (p^ from that predicted for rigid rods. However, we believe these 
theories to be accurate, specifically because simulations by Poison and Frenkel of hard semi- 
flexible rods demonstrate that flexibility increases the N-S transition concentration above 

n 

that predicted for rigid rods j3[. 

A second possible explanation for the discrepancy between our experimental values of ^f^^^. 
and the predicted values for the phase transition of semifiexible hard rods is that the elec- 
trostatic interactions between the rods are indeed significant at high ionic strength, and that 
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they are different from the predicted electrostatic interactions. As our measurements suggest 
that both the second virial and avoidance approximations for the electrostatic interactions 
describe a N-S transition which differs significantly from the experimentally observed be- 
havior, we believe this is possible. Further evidence to suggest that electrostatics still plays 
a role at high ionic strength is visible in Fig. El which shows that the surface charge of the 
rods can still influence the phase behavior, even when the phase transition has become ionic 
strength independent. 

A third possibility is that coupling electrostatics and flexibility produces an inter-rod 
repulsion which is a complex combination of flexible-hard rod and charged-rigid rod inter- 
actions. It has been shown for concentrated suspensions of DNA, that fluctuations due to 
the flexibility of the DNA actually enhance inter-rod repulsions in an exponential manner 
|85l | . The consequence is that for a given osmotic pressure exerted by a concentrated DNA 
suspension, the volume fraction of DNA is much lower than predicted by Poisson-Boltzman 
electrostatics alone. This hypothesis is consistent with our measurements of a (p^.^^ that is 
much lower than predicted for both rigid and flexible rods. 

As our results are quite unexpected with respect to current theoretical predictions for the 
nematic - smectic transition of rod suspensions, we have speculated as to possible explana- 
tions for our observations. To fully understand the phase behavior of charged, flexible rods 
further computational and theoretical work is clearly needed. 

VI. CONCLUSION 

We have examined the nematic-smectic phase diagram for charged, semiflexible virus 
rods as a function of length, surface charge and ionic strength. We observed that in the 
semiflexible-rod limit the N-S phase boundary is independent of rod flexibility, as predicted 
theoretically. However, by studying the ionic strength dependance of this transition we ob- 
served that renormalizing the measured phase transition volume fraction, 0'^, by Onsager's 
effective diameter, Des, does not produce an ionic-strength independent phase transition 
concentration. Therefore the second virial approximation cannot be used to map the mea- 
sured nematic - smectic phase transition of charged, flexible rods onto hard, flexible rod 
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theories. 

At high ionic strength, we found that the concentration of the N-S phase boundary is 
independent of ionic strength. Kramer and Herzfeld's avoidance diameter theory [7] quahta- 
tively reproduces the observed ionic strength independent N-S phase behavior, but predicts 
that the ionic strength independent phase boundary is equal to the predicted hard-rod phase 
boundary. Our experimental results, however, are three times lower than the hard-particle 
phase boundary predicted for semiflexible rods, 0^ J3|, |9|. Clearly, more theoretical work 
is needed to understand the nematic-smectic phase transition of charged, flexible rods, in 
order to reconcile the differences we observe with charged, flexible viruses and that reported 
in simulations of hard, flexible rods. 

Finally, significant differences were measured between M13 and fd nematic-smectic phase 
transition concentrations, even when they shared the same average surface charge. These 
results indicate that the electrostatic interactions between these rods are more complicated 
than can be accounted for by calculating the interparticle potential assuming a uniform 
renormalized surface charge. We hypothesize that the electrostatic interactions between 
rods is influenced by the configuration of the charged amino acids on the viral surface. Ex- 
perimental tests of this hypothesis could be made by measuring M13 and fd equations of 
state (pressure vs density), and thus the particle-particle interactions, as a function of solu- 
tion salt and pH, as in techniques developed for DNA i35'. Computationally, this hy pot hesis 
can be tested by calculating the pair potential between rods with discrete charges J3l|. 
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